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Abstract 

pg ■ Even though the "free" neutrons in the inner crust of a neutron star are superfluid, they are 

< | ....... 

still strongly coupled to nuclei due to non-dissipative entrainment effects. These effects have been 

systematically studied in all regions of the inner crust in the framework of the band theory of solids. 

Using concepts from solid-state physics, it is shown that the density of conduction neutrons, i.e. 

neutrons that are effectively "free" , can be much smaller than the density of unbound neutrons 

(by an order of magnitude in some layers) due to Bragg scattering. These results suggest that a 

revision of the interpretation of various observable neutron-star phenomena may be necessary. 
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Born from the catastrophic gravitational core collapse of massive stars with M ^ 8M 
at the end point of their evolution during type II supernova explosions, neutron stars are 
among the most compact objects in the Universe. The surface of a neutron star, which is 
composed mainly of iron, is generally obscured by a thin atmosphere. A few meters below 
the surface at densities above ~ 10 4 g-cm -3 , matter is so compressed that atoms, which 
are arranged in a regular Coulomb lattice, are fully ionized and coexist with a degenerate 
electron gas. Deeper in the star, nuclei become more and more neutron-rich ill . At densities 
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which is expected to be superfluid at low enough temperatures [3|]. The crust of a neutron 
star extends up to about ~ 10 14 g-cm -3 (i.e. half the density found inside heavy nuclei), at 
which point nuclei dissolve into a uniform mixture of electrons and nucleons. 

Various astrophysical neutron-star phenomena like pulsar glitches or quasiperiodic os- 
cillations in soft-gamma repeaters are expected to be related to the dynamical properties 
of the neutron superfluid in the inner crust (see e.g. [4j for a recent review). In model- 
ing these events, the neutron superfluid is generally assumed to be flowing freely through 
the crust because of the absence of viscous drag. This assumption is actually unrealistic. 
Even though the neutron superfluid can flow without friction, it can still be entrained by 



the crustal nuclei 



01? 



5K7J]. Unlike viscous drag, this entrainment effect is non-dissipative and 
therefore it persists in the superfluid phase. Entrainment effects in neutron-star crusts has 
been already estimated, but only in the shallow layers of the inner crust around the neutron 
drip density ~ 4.7x 10 11 g-cm -3 and near the crust bottom at densities above 5 x 10 13 g-cm -3 



using different crust models [8H10J . These results suggested that there might exist regions at 
intermediate densities where the neutron propagation could be completely suppressed owing 
to the presence of a band gap in the energy spectrum of unbound neutrons. Such crustal 
layers would resist the flow of a neutron current in the same way as an ordinary insulator 
resists the flow of an electric current. 

In this paper, entrainment effects are systematically studied in all regions of the inner 
crust of a cold non-accreting neutron star using a unified treatment based on the nuclear 



energy density functional theory (EDF), see e.g. HI] for a review. In particuliar, the 
static long-wavelength neutron current-current correlation function is computed for various 
densities, ranging from neutron drip to the crust-core transition. This function, which can 
be expressed in terms of an effective neutron mass or less ambiguously in terms of a density 



of conduction neutrons, is a necessary microscopic ingredient for a realistic hydrodynamical 
description of the neutron superfluid in the crust [5|, |6j . Both the equilibrium composition 
of the crust and the current- current correlation function have been calculated consistently 
using the same nuclear energy density functional. 

I. NUCLEAR BAND THEORY OF NEUTRON-STAR CRUSTS 



Since the seminal work of Negele&Vautherin [121 ]. the inner crust of a neutron star has 
been generally studied in the framework of the Wigner-Seitz approximation 13( according to 
which the crust is divided into a set of independent spheres centered around each lattice site. 
Each cell can thus be seen as an isolated giant "nucleus" for which the usual methods from 
nuclear physics can be applied. Even though this approach has been fruitful for calculating 
ground-state properties (at least at not too high densities, as discussed in Refs. [14J, [l5|), it is 
inappropriate for studying the low-energy dynamics of "free" neutrons which are delocalized 
over the whole crust like "free" electrons in ordinary metals. Indeed, the interactions of 
unbound neutrons with the crystal lattice are highly non-local thus leading to long-range 
correlations which cannot be properly taken into account in the Wigner-Seitz approximation. 

In this work, the neutron conduction in the neutron-star crust will be studied using 



the band theory of solids 16 1. Although this theory has been very successfully employed 



in various systems in optics, acoustics and condensed matter physics, its application to 



neutron-star crusts is rather recent |8Hl0j. The band theory relies on the assumption that 
the solid can be treated as a perfect crystal. Although the crust of a real neutron star might 
not be a perfect crystal (see e.g. Section 3.4 of Ref. [4j and references therein), this is still 
a reasonable approximation for cold non-accreting neutron stars. In the rest frame of the 
crust, both bound and unbound neutrons are supposed to be described by static periodic 
mean-fields. In the model we employ here, these fields are generated self-consistently using 
the nuclear energy density functional theory with a semi-local functional of the Skyrme 
type [ll|. This functional is of the form 

E = -Ekin + -^Coul + -^Sky , (1) 

where -Ekin is the kinetic energy of the sample volume, -Ecoui is the Coulomb energy (using 



the Kohn-Sham approximation for the exchange part [17|) and E Sky = f d 3 r £ S]iy (r) is the 



nuclear energy. Ignoring pairing, which represents a small correction to the total energy, 
.Esky is a functional of the following local densities and currents with q = n,p for neutron, 
proton respectively: (i) the number density 

««(»•) = Yl n g( r ' CT ' r ' a ) ' ( 2 ) 

<r=±l 

(ii) the kinetic density 

r «( r ) = E [ d 3 r'S{r-r')V-Vn q {r,a-y,a) , (3) 

and (iii) the spin current vector density 

•W = "~ E /dV5(r-r')(V-V')xn g (r,a;r',cr')(a'|a|a) , (4) 

where n q (r,a;r f ,a') is the density matrix in coordinate space (denoting the spin states by 
a, a' = 1,-1 for spin up, spin down respectively). Introducing the isospin index t — 0, 1 for 
isoscalar and isovector quantities respectiveljo, the Skyrme functional Esk y can be expressed 
as 

£ Sky = J] C>* + C t An n t An 4 + C^ t r t + C t VJ n t V • J t + \c(j t 2 , (5) 

i=0,l 

The C-coefficients are determined by fitting experimental nuclear data and/or properties of 
infinite homogeneous nuclear matter as obtained from many-body calculations. The coeffi- 
cients C" are generally not constant but depend on the (isoscalar) density n. Historically 
the Skyrme functional was obtained from the Hartree-Fock approximation using zero-range 
effective nucleon-nucleon interactions of the kind 

v itj = t (l + x P ff )6(rij) + -ii(l + xiP ff )^-2 [p% 5(r tJ ) + 6^)$] 
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+ -2 H^o (£< + ^i ) • P« X ^ (Tij ) Pij , (6) 

where r^ = r j — r.,-, r = (r, + fj)/2, pij = —ih(Vi — Vj)/2 is the relative momentum, P CT is 
the two-body spin-exchange operator. Minimizing the total energy E for fixed numbers of 
neutrons and protons leads to the set of self-consistent equations 

Mr)^(0 = eiV2(r), (7) 



1 Isoscalar quantities, also written without any subscript, are sums over neutrons and protons (e.g. no = 
n = n n + rip) while isovector quantities are differences between neutrons and protons (e.g. n\ — n n — n p ). 



in which the single-particle (s.p.) Hamiltonians h q (r) can be expressed as 

hq (r) = -V ■ B q (r)V + U q (r) - iW q (r) ■ V x a , (8) 



with the various fields defined by 
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These fields depend on the local densities and currents, Eqs. ()2]), ([3]) and (j4]), which in turn 
can be expressed in terms of occupied s.p. wavefunctions f^l(r) (see e.g. Ref. H[). The 
periodicity of the crystal lattice means that 

B q (r + t) = B q (r) 1 (10) 

U q (r+£) = U q (r), (11) 

W q {r+l)=W q {r), (12) 

for any lattice translation vector £. The boundary conditions to be used in Eqs. ©-© are 



imposed by the Floquet-Bloch theorem 16] 



^l(r+£)=ex P (ik-e)^ ) k (r), (13) 

where k is the Bloch wave vector and a is the band index. 

Because protons are tightly bound to nuclei, proton band structure effects are very small 
(i.e. the proton s.p. energies are essentially independent of k) and will not be discussed 
here. In the following, the neutron energy bands and neutron wave functions will thus be 
denoted simply as e a k and (p a k(f) respectively. 

II. CONDUCTION NEUTRONS 

The static long-wavelength neutron current-current correlation function determines the 
mass current of the neutron liquid induced by a change of crystal momentum, the lattice 
being fixed. Neutron pairing, which gives rise to neutron su per fluidity, is expected to have 
a minor impact on the current-current correlation function 18j and will thus be neglected 
for simplicity. In the ground state, all neutron s.p. states lying below the Fermi level are 
occupied so that the neutron liquid is at rest in the crust frame. Considering a small shift 



5k in fe-space of the Fermi surface (FS), each s.p. state will then carry on average a net 
crystal momentum p n = h5k leading to a net mass current j n . Neglecting back flow effects 
(i.e. assuming that the s.p. energies remain unaffected), the current is given to first order 
in 5k by 

d 3 k 



Jn = m. 
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5h ak V k e a k = n c n p n 



(14) 



where m n is the neutron mass, 5n a k = n a k — n^ k denotes the change in the distribution 
function {n° ak being the ground-state distribution) while the density rf n is defined in terms 
of the trace of the effective mass tensor introduced in solid-state physics [16J 
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Incidentally the neutron effective mass tensor (15) has been also introduced for the study 
of neutron diffraction in ordinary crystals 19|, |20|]. Since the ground-state distribution is 
simply given by fiP ak = H(ep — s a k), where H{x) is the Heaviside unit-step distribution, 
the integral in Eq.( TT6l) has to be taken over the Fermi volume. The Fermi energy e-p is the 
Lagrange multiplier introduced during the minimization of the total energy ([T]) in order to 



ensure the conservation of the neutron number 



n, 



E 



d 3 A; 

(2^ 



r,0 



-n 



ak 



(17) 



The density n c n can be interpreted as the density of conduction neutrons by analogy with 
conduction electrons in ordinary metals. These conduction neutrons are not entrained by 
nuclei and can thus be considered as being effectively "free" . Equation (|14p implies that in 
an arbitrary frame where the crust moves with a velocity v p , the neutron mass current will 
no longer be aligned with the neutron momentum, but will be given by 



jn = n c n p n + {n n - n c n )v p . 



(18) 



The quantity n n — rf n can be interpreted as the density of neutrons that are effectively bound 
to nuclei. 



Alternatively, the neutron conduction can be expressed in terms of an effective mass m* 
by writing the neutron momentum in the crust frame 

p n = m>„ , (19) 

where v n is the average velocity of free neutrons defined by 

j n = n { n m n v n . (20) 

Comparing Eqs. ffl9l and (120]) with (fT4"j) leads to the following expression for the effective 

mass 

f 
n 

< = m n ^ . (21) 

III. NUMERICAL RESULTS 

In principle, the ground-state structure of the neutron-star crust at a given average nu- 
cleon density could be determined by solving self-consistently the EDF equations ©,([8]) 
and (Q with Bloch boundary conditions ( TT3|) under the constraint of beta equilibrium. Con- 
sidering that the equilibrium structure of the crust is a body centered cubic crystal, these 
calculations should be repeated for different lattice spacings until the lowest total energy (CQ) 
is found. Such calculations would be computionally extremely expensive. For this reason, 
we have taken the composition of the crust, as found in Ref. [2|. These calculations were 
based on the fourth-order Extended Thomas-Fermi method with proton shell effects added 
via the Strutinsky-Integral theorem. Neutron shell effects were neglected since they were 



shown to be much smaller than proton shell effects 2l| . This so-called ETFSI method is 



actually a high-speed approximation to the self-consistent EDF equations (Ej),® and (|9]). 



p 



The functional BSkl4 22J used in Ref. |2( was not only fitted to essentially all the available 
experimental atomic mass data with a root mean square deviation of 0.73 MeV, but was also 
constrained to reproduce the neutron- matter equation of state of Ref. [231 ]. obtained from 
many-body calculations using realistic two- and three-body nucleon-nucleon interactions. 
As a matter of fact, this equation of state is in good agreement with more recent calcula- 



tions 



2J-I26 



261 ] in the density domain relevant to neutron-star crusts. For all these reasons, 
the Skyrme interaction BSkl4 is particularly well-suited for studying neutron-star crusts. 
The neutron and proton density distributions are shown in Fig. [1] for a few crustal layers. 
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FIG. 1. Neutron (solid line) and proton (dashed line) density profiles inside the Wigner-Seitz 
cell for different baryon densities n (in fm -3 ), as obtained with the ETFSI method [2]. Note the 
formation of "bubbles" at n = 0.08 fm -3 : the nucleon densities are slightly larger at the cell edge 
than at the cell center. 

The neutron energy bands e a k have been calculated by solving Eqs. (CJ),® and ([9]) with 
Bloch boundary conditions (TT31) using the self-consistent fields obtained in Ref.[2j. The 
spin-orbit potential W n (r), which is proportional to Vn n (r) and Vn p (r), is much smaller 
in neutron-star crusts than that in isolated nuclei and have therefore been neglected. The 
neutron band structure has been computed by expanding the neutron s.p. wave functions 
into plane waves 



Lp k (r) = exp(ifc ■ r) ^ <Pk(G) exp(iG ■ r) 



(22) 



G 



where G are reciprocal lattice vectors. By definition, exp(iG • t) = 1 so that the Bloch 
boundary conditions f fl~3|) are automatically satisfied. In this way, the EDF equation for the 
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FIG. 2. Left panel: neutron band structure in the inner crust of a neutron star at the average 
baryon density n = 0.0003 fm~ 3 along high-symmetry lines in the first Brillouin zone (only unbound 
states are shown) . The arrow indicates the position of the neutron Fermi energy. Right panel: band 
structure of a uniform neutron gas at density n^ (reduced zone scheme) . For comparison with the 
left panel, all bands have been slightly shifted. 



neutrons reduces to the matrix eigenvalue problem 

(* + G a ) ■ (k + G„)B n (Gg -G a ) + V n (G f -G a ) 



E 



¥k{Gp) = e ak ip k {G a ) (23) 



with 



B n (G) 



d 3 r J B n (r)exp(-iG-r), U n (G) 



1 



d 3 rl/" n (r)exp(-iG-r). (24) 



Kell J Kell 

These integrals, taken over any primitive cell of volume V ce n, can be efficiently calculated 
using fast Fourier transforms. 



In the shallow layers of the inner crust, the linearised augmented plane wave method 10 1 
would have been computionally much faster. Unfortunately this method can be reliably 
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FIG. 3. Left panel: neutron band structure in the inner crust of a neutron star at the average 
baryon density n = 0.03 fm -3 along high-symmetry lines in the first Brillouin zone (only unbound 
states are shown). The arrow indicates the position of the neutron Fermi energy. Right panel: 
band structure of a uniform neutron gas at density n n (reduced zone scheme). For comparison 
with the left panel, all bands have been slightly shifted. 



applied only in the vicinity of the neutron drip transition where only a few neutron bands 
(in the continuum) are filled. For this reason, the plane wave method has been used in all 
regions of the inner crust. 

The number of calculated neutron bands varies from a hundred in the shallowest layers of 
the inner crust at n = 0.0003 up to 850 at n = 0.03 fm~ 3 . A few band structures are shown 
in Figs I2"f3l and HI We have also shown the band structure in the reduced zone scheme 16| of 
a uniform neutron gas of density n n (empty lattice limit), whose energies are simply given 
by 



e ak = B n (n f n ){k + G a ) 2 + U n (n l n ). 
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FIG. 4. Left panel: neutron band structure in the inner crust of a neutron star at the average 
baryon density n = 0.08 fm -3 along high-symmetry lines in the first Brillouin zone (only unbound 
states are shown). The arrow indicates the position of the neutron Fermi energy. Right panel: 
band structure of a uniform neutron gas at density r^ n (reduced zone scheme). For comparison 
with the left panel, all bands have been slightly shifted. 

As can be seen in Fig [2J the energy spectrum of unbound neutron states near the neutron 
drip point very closely ressembles that of a uniform neutron gas, even though the periodic 



mean-field potential drops by ~ 50 MeV inside nuc 



of the Phillips-Kleinman cancellation theorem 27|. The orthogonalization of the unbound 



ei. This striking result is the consequence 



states to the bound states leads to an effective repulsive non-local and energy-dependent 
potential which partially cancels the strongly attractive mean-field potential of nuclei. To a 
large extent, the same kind of cancellations occurs in the densest regions of the crust as can 
be inferred from Fig HJ On the contrary, the neutron band structures in the intermediate 
layers of the inner crust differ significantly from that of a uniform neutron gas, as illustrated 
in Fig [3j In particular, the energy spectrum is much denser thus revealing that on average 



11 



the energies e a k have a much weaker k- dependence (almost flat bands) than that given by 
Eq. fl25|) . As a result, the conduction neutron density n c n given by Eq. f|T6|) is expected to 
be much smaller than the density n n of unbound neutrons at densities n ~ 0.03 fm -3 . 

The computation of the conduction neutron density requires the evaluation of the effective 
mass tensor (1151) for all occupied bands. In fact, it follows from Eq. ( II 6p and the periodicity 
of the s.p. energies that completely filled bands do not contribute to the current. This can 
be easily seen by simply rewriting Eq. (fT6|) as an integral over the Fermi surface using the 
Green-Ostrogradsky's theorem 

m r , 



nl 



^E/jV^d^. ( 26 ) 



24tt 3 ^ 

Even though this expression is fully equivalent to Eq. (TT6|) . it is computationally much more 
convenient since only the evaluation of the first derivative of e a k is needed. In addition, these 



derivatives can be easily calculated analytically using the Hellmann-Feynman theorem 28] 

o± = j2 ^ G «y b ^°p - G ^ 2k% + G » + G ti &( G /o ( 2? ) 

G a ,Gp 

with the wavefunctions normalized as 

^|^(G,)| 2 = 1. (28) 



For each average density n, the neutron Fermi energy Ep has been determined solving Eq. ([1 



using the mean- value point method 29J. The Fermi surface integral in Eq. (fT8l) has been 
evaluated with the Gilat-Raubenheimer method [3 Of using up to 1360 points in the irre- 
ducible domain (i.e. 65280 points in the first Brillouin zone) in order to ensure a precision 
of a few percent. Results are summarized in Table [B As expected from the band structures, 
the flow of neutrons is almost unaffected by nuclei in the peripheral regions of the inner 
crust. On the contrary, the neutron conduction is found to be almost completely suppressed 
at densities n ~ 0.02 — 0.03 fm~ 3 . Whereas more than 90% of neutrons are unbound in these 
layers, less than 10% of them are actually conducting leading to a huge enhancement of the 
neutron effective mass m* ~ 13.6m n . Incidentally, this result is in close agreement with 
the effective mass m s n tar ~ 15.4m n obtained in a previous work [9] using a different crust 
model thus suggesting that such strong entrainment effects are generic. However, further 
work remains to be done exploring the dependence of m* on the nuclear energy density 
functional. 
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n (fm 3 ) 


Z 


A 


ni/n n (%) 


</< (%) 


m*/m n 


0.0003 


50 


200 


20.0 


82.6 


1.21 


0.001 


50 


460 


68.6 


27.3 


3.66 


0.005 


50 


1140 


86.4 


17.5 


5.71 


0.01 


40 


1215 


88.9 


15.5 


6.45 


0.02 


40 


1485 


90.3 


7.37 


13.6 


0.03 


40 


1590 


91.4 


7.33 


13.6 


0.04 


40 


1610 


88.8 


10.6 


9.43 


0.05 


20 


800 


91.4 


30.0 


3.33 


0.06 


20 


780 


91.5 


45.9 


2.18 


0.07 


20 


714 


92.0 


64.6 


1.55 


0.08 


20 


665 


104 


64.8 


1.54 



TABLE I. Composition of the inner crust of cold non-accreting neutron stars as obtained from 

Q 

Ref. [2j. Z and A are respectively the average number of protons and the total average number of 
nucleons inside the Wigner-Seitz cell. n n is the average neutron density, n n is the density of free 
neutrons as defined by the quantity psn in Ref. |2j , n c n is the density of conduction neutrons and 
m^ the neutron effective mass. Note that in the densest layer, n n > n n due to the formation of 
bubbles as indicated in Fig. [Q 



IV. MICROSCOPIC ORIGIN OF ENTRAINMENT 



The large discrepancy between the density of unbound neutrons and the density of con- 
ducting neutrons is somehow counterintuitive. Indeed in ordinary metals, the electrons that 
are tighly bound inside the individual atoms constituting the solid have their wavefunction 
vanishing exponentially outside atoms and are therefore not much affected by the Bloch 
boundary conditions. As a consequence, their energy bands in fc-space are essentially flat so 
that Vfc£ Q A: — hence yielding a negligible contribution to the current. The non-trivial elec- 
tron band structure arises from the most loosely bound "valence" electrons in the isolated 
atoms which become delocalized in a metal and which can be generally identified with the 
conduction electrons (still, the density of valence electrons is not exactly equal to the density 
of conduction electrons). On the contrary, the neutron-saturated "nuclei" found in the inner 

13 



crust of a neutron star only exist because of the Pauli blocking effect from the surrounding 
neutron liquid but would decay immediately in vacuum. For the reasons mentioned above, 
neutrons bound inside nuclei do not contribute to the current. Since unbound neutrons are 
delocalized, one might naively expect that they are all conducting. Indeed, ignoring the 
crystal lattice and treating the unbound neutrons as a uniform gas of density n n , it follows 
immediately from Eq. (fl6|) or Eq. (|26|) that n c n = n n . However it should be emphasized 
that the density of conduction neutrons is fundamentally different from the density of un- 
bound neutrons: the former characterizes the dynamics of the neutron liquid while the latter 
is a ground-state property. These two densities are generally not equal because unbound 
neutrons can be scattered by the crystal according to Bragg's law. 

The effects of Bragg scattering are embeddied in the effective mass tensor ffT5|) appearing 
in the definition ( Tl6|) of the conduction neutron density. The components of this tensor need 
not be positive and can actually be negative for wave vectors k such that Bragg reflection 
is allowed. In particular, an unbound neutron will be reflected whenever its Bloch wave 
vector k lies on a Bragg plane, i.e. k satisfies the von Laue condition Ik ■ G = G 2 for any 
reciprocal lattice vector G. As a result, the neutron current will be suppressed however small 
the periodic mean- field potential is. Using perturbation theory and setting B n (r) = B n for 
simplicity, it can be shown that in the vicinity of a Bragg plane the s.p. energies of free 



neutrons are split into two bands [16] 



4 ~ l -B n {k 2 + (k- Gf) ± ^\&JP - (* " G) 2 ) 2 + Pn(GW ■ (29) 

It is easily seen that the effect of the periodic potential is to flatten the bands around Bragg 
planes by introducing a gap of magnitude 2\U n (G)\. As a consequence, the conduction 
neutron density will be reduced whenever the Fermi energy lies on the top (bottom) of the 
lower (upper) band. The neutron conduction thus depends on two factors: i) the amount of 
flattening of the bands which in turn is governed by the periodic potential, ii) the position 
of the Fermi level which is determined by the total neutron density. The more the Fermi 
sphere intersects Bragg planes, the larger will be the effect of Bragg scattering on the neutron 
conductivity. The number of intersections depends on the ratio between the Fermi volume 
and the volume of the first Brillouin zone. According to Luttinger's theorem [31[, the Fermi 
volume associated with unbound states is given by Vp = (2iT) 3 n t n . By definition, the volume 
of the first Brillouin zone is equal to Vbz — (27r) 3 /Vccii- Therefore the ratio Vf/^bz = nlYceii 

14 



is simply equal to the average number of unbound neutrons inside the Wigner-Seitz cell 
(or equivalently to the average number of unbound neutrons per nucleus). Basically, this 
number is the lowest at the neutron drip point, peaks at about 1417 at density n = 0.03 
fm -3 and decreases at higher densities. As expected, the resistance of the crust to the flow 
of a neutron current follows a similar behavior (see Table H]). 

V. CONCLUSIONS 

Despite the absence of viscous drag, the neutron superfluid permeating the inner crust of 
a neutron star is still strongly coupled to nuclei due to non-dissipative entrainment effects. 
These effects have been systematically studied in all regions of the inner crust of a cold 
non-accreting neutron star in the framework of the band theory of solids using an effective 
Skyrme nuclear energy density functional which not only yields an excellent fit to essentially 
all experimental atomic masses but was also constrained to reproduce a realistic neutron- 
matter equation of state. 

Some regions of the inner crust are found to strongly resist the flow of a neutron current in 
the same way as an ordinary insulator resists the flow of an electric current. As a consequence 
the density of conduction neutrons, i.e. neutrons not entrained by nuclei, turns out to be 
much smaller than the density of unbound neutrons thus leading to a huge enhancement 
of the neutron effective mass. These results suggest that a revision of the interpretation of 



many observable astrophysica 



phenomena like pulsar glitches |32|, quasiperiodic oscillations 
in Soft-Gamma Repeaters J33j and the cooling of neutron stars 34J, may be necessary. 

Incidentally, the kind of entrainment effects discussed here are also expected to be ob- 
served in laboratory superfluid systems. Indeed, qualitatively similar results have been 
predicted for a unitary Fermi gas in a ID optical lattice J35(. Since a dilute neutron gas is 
approximately in the unitary regime [36|, cold atoms experiments could shed light on the 
properties of superfluid neutrons in neutron-star crusts. 
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